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ABSTRACT  (ONGERUBRICEERD) 

Atmo^henc  structures  can  be  measured  widi  incoherent  (q>tical  ladars  Qidars).  Because  these  structures 
drift  with  the  wind,  they  can  serve  as  a  tracer  for  remote  sensing  of  the  wind  vector.  For  this  purpose,  a  dual 
mraostatic  sr-mwiing  lidar  system  is  available  to  measure  die  atmosphere  simultaneously  in  two  different 
directions  over  a  maximum  range  of  about  1  km.  The  transit  time  of  identified  patterns  between  two  sensing 
points  in  die  horizontal  plane  provides  in  combination  widi  die  geometry  of  the  lidar,  sufficimt  information 
to  derive  die  horizontal  wind  vector.  The  method  is  based  <«  cross-correlaticHi  techniques.  To  determine  die 
spatial  wind  vector  as  a  function  of  aldtude  it  is  sufBcient  to  measure  in  duee  different  upward  directions. 
This  can  be  realized  with  a  trple  lidar  or  widi  a  single  lidar  measuring  ctmsecutivety  in  duee  different 
directions  and  using  an  equivalent  but  more  extended  inversion  mediod. 

This  rqxvt  describes  die  main  activities  carried  out  for  die  project  DWARSWIND'  in  the  period  January 
1991  -  Mardi  1992.  After  a  summary  of  a  short  literature  study,  some  dieoretical  aqiects  ate  described  such 
as:  die  vector  rqnesaitation  of  fhe  dual  hdar  in  Cartesian  cotndinates,  a  mediod  to  derive  dte  wind  vector, 
die  characteristic  life  time  and  die  diaracteristic  size  of  the  structures  frun  a  set  of  lidar  measurements.  A 
setectkm  of  experimental  results  are  presented.  For  exanqde,  die  cakulatioii  of  die  honzcmtal  wind  vector  at 
an  ahitude  of  15  m  and  die  wind  vector  iq>  to  an  aUtnde  of  1000  m. 
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SAMENVATTING  (ONGERUBRICEERD) 

Met  behu^  van  een  incoherente  lidar  (optische  radar)  kunnen  atmosferische  strukturm  gedetecteerd  wordea. 
Omdat  deze  strukturen  met  de  wind  meedrijven  Imnnen  ze  gebraikt  wordeo  als  hu^nniddel  voor  bet  op 
afstand  bepalm  van  de  windvector.  Eoi  monostatiadi  tweevoudig  lidarsystem  is  beschikbaar  om  gelijktijdig 
in  twee  verschillende  richtingen  te  kunnen  metoi  over  een  beieik  van  maxhnaal  1  km.  Uit  de  loopdjd  van 
herkoobate  aerosol-strukturen  tusaen  twee  meeqxinten  in  bt;t  hcnizontale  vlak  en  de  geometiie  van  de  lidar, 
kan  met  bdiulp  van  kruiscoirelatie  de  windvector  wtmlen  bqpaald  in  bet  vlak  waarin  de  metingen  wotd«i 
uitgevoeid.  Om  de  windvector  als  fiinctie  van  de  boogte  te  bqMlen  is  bet  voldoende  om  in  drie  verschillende 
richtingen  sdinin  omboog  te  meten.  Dit  kan  geiealiseerd  woiden  met  een  drievoiidige  lidar  of  met  eoi 
enkelvoudige  lidar  die  achteieenvolgens  in  drie  verschilkode  richtingen  meet 

Dit  nqpport  beschrijf)  de  boofdactiviteiten  zoals  uitgevoeni  in  bet  kader  van  bet  project  'DWARSl^lND'  in 
de  periode  januari  1991  -  maart  1992.  Na  een  bdmqtte  litentuor  stndie  wortei  enige  theoretische  aq)ecten 
besdneven  zoals  een  vectoriSle  voorslelling  van  de  lidar  in  bet  Cartesuanse  assenstelsd  en  een  m^ode  om 
uit  een  set  lidar  meringen  de  windsnelheid  te  berekenen.  Tevens  wordt  aangegeven  hoe  van  de  atmosferische 
strukturen  de  karakteristkke  grootte  en  de  karakteristieke  levensduur  bq>aald  kunnen  wordoL  Tenslotte 
warden  enige  resultaten  gqireseineeni  die  met  de  lidar  zijn  vertaegen,  zoals  de  horizontale  windvector  op 
een  boogie  van  ca.  ISmendemimtelijkewindvector  tot  <y  een  boogte  van  maximaal  1000  m. 
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1  INTRODUCTION 

Wind  is  an  in^witant  meteorological  parameter,  which  is  generally  measured,  e.g.  by  in  situ 
mechanical  and/or  acoustical  sensors  (ultrasonic  anemometers).  In  those  situations  where  the 
wind  vector  as  a  function  of  altitude  is  required,  for  in  stance  for  die  prediction  of  ballistic 
trajectories,  balloon  measurement  can  provide  the  required  infonnatioo.  Cmrently,  it  is  also 
possible  to  use  remote  sensing  like  Dt^kr  radar  (e.g.  Vaisala  Model  400)  or  lidar  to 

measure  the  wind  vector.  The  advantage  of  lidar  is  diat,  in  die  future,  mote  tasks  can  be  cmnbined 
within  one  qitical  system.  For  instance,  for  measuring  range,  visibility,  vertical  structure  of  die 
atmospheric  exdncdon  and  wind  vector  as  well. 

It  has  been  shown,  among  odiers,  by  Zuev  et  al.  (1973),  Eloranta  et  al.  (1975),  Sasano  (1985), 
Eloranta  and  Schols  (1990)  and  Kunz  (1991),  diat  it  is  possible  to  measure  dynamic  atmos^eric 
structures  with  incrdieroit  optical  radars  (lidars).  If  it  is  assumed  that  these  structures  drift  with 
the  wind,  they  can  serve  as  tracers  for  wind  measurements.  The  atmosidieric  structures  are 
characterized  by  a  diffetenr  aerosol  concentration  and/or  by  a  different  refractive  index  with 
respect  to  its  surroundings.  The  sizes  of  dtese  structures  vary  from  millimetres,  as  proven  by 
Frehlich  (1988)  to  kilometres  as  shown  by  Hardy  (1969),  Konrad  (1970),  Zuev  et  al.  (1977), 
Sasano  (1985),  Hooper  and  Eloranta  (1986)  and  Kolev  et  al.  (1988).  Life  times  of  the  structures 
are  in  die  range  from  a  few  seconds  to  tens  of  seconds  and  more  as  measured,  e.g.,  by  Zuev  et  al. 
(1977)  and  Derr  et  al.  (1979).  Balin  (1990)  mentioned  diat  die  size  of  the  structures  depends  (m 
the  wind  speed  and  on  the  tenqierature.  It  has  also  been  shown  e.g.  by  Hardy  (1969),  Konrad 
(1970),  Batten  (1973),  Noonkester  (1976)  and  Brookner  (1977)  that  these  structures  can  be 
detected  radar. 

The  princ^le  to  derive  the  wind  vector  frtnn  a  set  of  incoherent  lidar  data  is  best  illustrated  in  a 
one-dimeosional  qpace.  In  diat  case,  die  structures  move  along  a  line  and  their  velocity,  dieir 
duracteristic  life  time  and  dieir  characteristic  size  can  be  inverted  from  the  time  bdiaviour  of  die 
structures  measured  at  two  different  locatioiis  (note  that  one  range  bin  of  a  lidar  wavdfonn  is 
equivalent  with  me  point-sensmr).  The  cross^orrelation  <ji  die  two  time  series  and  die  separation 
between  die  sensors  provides  a  direct  measure  for  die  parameters  mentioned,  fr  is  evident  dud  the 
quality  of  die  calculations  dqieiids  on  the  dynamics  (hfe-time)  of  the  drifting  structures. 

By  means  of  two  lidars,  ifriiich  operate  simultaneously  in  one  jdane,  die  moving  structures  can  be 
detected  subsequently  at  different  locations.  Tlie  disidacement  can  be  visualized  by  plotting  die 
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data  in  a  rima  versus  range  nuq>  (assuming  diat  die  structures  move  <nily  in  this  plane).  Tbe  ratio 
of  die  vectorial  displacement  of  structures  and  the  time  difference  of  detection,  provides 

a  direct  measure  for  the  wind  vector. 

Instead  of  two  line-sensors  (lidars),  one  can  also  use  diree  point-sensors  to  determine  die  wind 
vector  in  a  plane.  In  dial  case,  die  geometry  of  die  sensors  and  die  mutual  cross-ctxrelatitm  of  die 
measured  time-series,  can  also  provide  the  wind  vector  in  diat  plane.  Hie  diree-point  method  can 
also  be  siniulated  with  one  lidar  by  measuring  subsequently  in  diree  different  directions.  This 
method  has  been  described  by  e.g.  Eloranta  (1990)  and  Matvienko  (1990).  CLvnesha  (1981)  noted 
that  this  mediod  has  already  been  applied  in  1949  to  map  moving  ionoqdieric  structures  at  radio 
frequencies.  Contrary  to  what  is  generally  assumed,  die  diree  sanqile  points  need  not  necessarily 
be  in  a  horizontal  plane,  as  will  be  shown. 

After  a  successful  lidar  experiment,  the  calculatitm  of  the  wind  vector  can  be  considered  as  a 
con^lex  signal-analyses  problem  with  large  amounts  of  data.  From  the  literature  two  different 
variants  of  inversion  methods  are  distinguished: 

1.  direct  correlation  of  the  structures  measured  at  two  or  three  positions  as  described  by  Briggs 
(1950),  Phillips  (1955)  and  Briggs  (1968).  Ferdinandov  (1982  and  1984)  took  into  account 
die  atmospheric  attenuation. 

2.  indirectly  by  transforming  the  range  versus  time  data  sets  into  wavenumber  vmus  time 
using  the  Fourier  transform.  The  drift  of  the  patterns  can  now  be  inverted  from  die  phase 
shift  of  the  spectral  conqionents  in  consecutive  lidar  returns  (cross  spectral  jdiase  shift).  This 
method  has  been  applied  by  Stroga  (1980)  and  Hooper  (1986).  The  madiematical  basics 
have  also  been  described  by  Stull  (1989). 

This  report  describes  the  main  activities  performed  in  the  first  period  of  die  project 
HWARSWIND'.  Geometrical  and  dieoretical  aqiects  ate  discussed  with  reqiect  to  die  dual  lidar, 
with  enqdiasis  on  calculation  of  die  wind  vector.  Some  typical  eiqierimental  data  will  be 
pres«ited  and  elaborated  to  illustrate  die  data  format  and  the  potential  of  the  system  for  wind 
measurement.  The  properties  of  die  lidar  in  question  have  been  described  in  an  earlier  rqiort  by 
Kunz  (1990).  Secondary  activities,  like  {nogramming  of  a  hardware  fari  Fburier  transformer  and 
an  analysis  of  die  noise  properties  of  die  receiver,  have  been  described  respectively  in  Kunz 
(1992a)  and  Kunz  (1992b). 
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2  VECim  REreESENTATION  OF  THE  DUAL  UDAR 

The  iniipose  of  this  section  is  to  define  the  vector  Rfiresentatioo  of  die  dual  lidar  in  Cartesian 
t  coordinates. 

A  siiiq>lified  line  model  of  the  monostatic  dual  lidar  ^stem,  with  die  two  beams  and  the  axes  of 
rotatim,  is  shown  in  Figm  2.1.  The  large  box  represents  die  main  lidar  and  the  mall  c»ie 
represents  the  auxiliary  lidar  which  has  been  mounted  on  top  of  the  main  lidar.  The  main  lidar  can 
rotate  around  an  horizontal  axis  £,  e.g.  for  adjustment  of  the  elevation  angle  during  vertical  scans. 
The  auxiliary  lidar  has  beat  mounted  on  die  mam  Udar  such  that  it  can  rotate  around  the  axis  p  for 
changing  the  angle  between  die  two  lidar  directimis.  Axis  p  is  perpendicular  to  the  upper  plane  of 
the  main  lidar.  The  two  lidar  beam  axes  g  and  pg  are  always  parallel  to  this  plane.  Tht  lual  lidar 
has  been  mounted  m  a  jack  which  can  rotate  around  the  vertical  axis  z  pointing  the  whole 
system  in  any  desired  azimuthal  direction.  By  scanning  both  in  elevation  and  in  arimuth,  the 
system  can  probe  the  whole  hemu^here. 


Rgnre  2.1:  Simplified  line  model  of  the  monostatic  dual  lidar  system  with  an  indication  of  die  basic 
vectors. 
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Rgine  'll  shows  die  rqiceseotatioD  of  die  dual  hdar  in  the  Cvtesum  cooniiiute  system.  The  x- 
axis  and  the  y-axis  define  the  hotizontal  plane.  The  z-axis  is  pqpendicular  to  this  ground  plane 
and  points  at  die  zenidi.  Note  diat  die  aTimiith  angle  in  die  Cartesian  oxHdinate  syston  is 
measured  counter  clockwiae  widi  respect  b>  die  x-axis  while  die  wind  directum  is  measured  clock 
wise  with  reflect  to  die  North.  The  elevation  angle  is  measured  widi  reqiect  to  die  ground  plane. 


Figure  2.2:  Vector  rqxesentation  of  most  major  axes  of  the  dual  lidar  system  and  d^nitioa  of  the 
reference  points. 

The  definition  of  the  vectors  and  die  angles  is  as  follows: 

A  s  unit  vector  in  die  measuring  direction  of  die  auxiliary  Udar 
£  ■>  unit  vector  of  die  elevatkm  axis,  parallei  widi  die  ground  plane 

m  *  unit  vectmr  in  die  measurii^  direction  of  die  main  lidar 

P  a  unit  vector  of  die  rotation  axis  of  die  auxiliaiy  lidar 
£  a  unit  vector  in  die  direction  of  die  wind 

A  8  unit  vectOT  in  the  diiection  of  the  Nordi 

y  8  unit  vector  in  die  direction  of  die  West 

2  a  unit  vector  in  the  vertical  direction 
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s  aziinudi  angle  of  tbe  main  lidar 
0^  s  etevaticm  angle  of  the  main  lidar 
0^  s  azimudi  angle  of  die  auxiliary  lidar 
0^  B  elevation  angle  of  die  auxiliary  lidar 
6  s  angle  between  main  lidar  and  auxiliary  lidar 


The  measuring  direction  of  die  main  lidar .  m.  is  given  by: 


m 


cos(0ea,)  •  COS(0_) 
cos(0_)  •  8in(0^ 

SUl(0eB,) 


(2.1) 


Because  m  is  expressed  as  a  unit  vector,  the  z  component  is  equal  to  the  sine  of  the  elevation 
angle  whereas  the  tangent  of  the  azimuth  angte  is  determined  by  the  quotient  of  the  x  and  the  y 
coiiqxment. 


The  axis  of  rotation  for  the  elevation,  £,  is  perpendicular  to  the  direction  of  the  main  beam,  qi,  and 
parallel  to  the  ground  plane.  Thus  one  can  write: 

£.ia  =  0  (2.2) 


The  coordinates  of  the  unit  vectOT  £  are  thoefore: 


sin(0«„) 
£  «  -cos(0^) 

0 


(2.3) 


The  measuring  direction  of  die  aiiTiliary  lid^  is  determined  the  measuring  direction,  m,  of 
the  main  lidar  and  the  angle,  6,  betweoi  die  two  lidar  systems.  Because  £  and  m  form  a 
rectangular  coordinate  system,  the  vector  a  can  be  expressed  as : 

a  =  m  •  cos(0)  +  £  •  sin(6)  (2.4) 


The  actual  direcdtHi  of  the  auxiliary  system  can  now  be  found  by  substituting  (2.1)  and  (2.3)  in 
(2.4).  This  results  in: 


a 


CO8(0)  •  co8(0g^  •  CO8(0^)  +  sin(0)  •  8in(0^ 
CO8(0)  •  cos(0  •  8in(0^  -  8in(0)  •  co8(0^ 
cos(0)'sia(0^) 


(2.5) 
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AddituHi  of  the  squared  CMupoamts  of  (2  J)  shows  that  g  is  indeed  a  unit  vectw.  As  a  result,  the  z 
conqKment  of  g  is  equal  to  the  sine  of  the  eievatitm  angle  of  the  auxiliary  system  whereas  (the 
tangent  of)  die  azimudi  angle  is  determined  by  die  quotient  of  the  x  and  die  y  conqionents.  Thus 
we  obtain  for  die  sine  of  the  elevation  angle  of  die  auxiliary  lidar. 

ain(0eii)  -  cos(0)  •  sin(0en])  (2.0) 

and  for  the  tangent  of  die  azunudi  angle  of  die  auxiliary  lidar  we  find: 

cos(0)  •  cos(0en,)  •  sin(0^)  -  sin(0)  •  cos(0^ 

tan(0O  =  - - - - — - SL.  (2.7) 

cos(0)  ■  cos(0eni)  ■  cos(0^  +  8in(0)  •  sin(0|^) 

The  axis  of  rotation  of  the  auxiliary  lidar,  |2,  perpendicular  to  the  tnain  lidar  is  found  from  the 
vector  product  of  g  and  q,  as  pointed  out  in  Hgure  2.2. 

E  =  ixm(*=fixiB)  (2.8) 


Substitution  of  (2.1)  and  (2.3)  in  (2.8)  provides  the  coordinates  for  the  vector  p: 


E  = 


sin(0^)  •  cos(0„_) 
sin(0^)  ■  sin(0^) 


(2.9) 


To  coiiqilete  Figure  2.2,  the  lines  ME  and  OB  are  calculated  with: 


m  +  Tr(g-Dl)  =  li£, 


(2.10) 


in  which  yand  p.  are  the  Imgth  of  respectively  ME  and  OE. 


Because  the  z-conqxment  of  g  is  zero,  the  value  of  y  can  be  found  by  siihatimting  the  z- 
conqxments  of  (2.1)  and  (2.5)  in  (2.10).  This  results  in: 


Y 


1 

1  -  CO8(0) 


(2.11) 


Substitution  of  (2.1),  (2.5)  and  (2.11)  in  (2.10)  {vovides  the  following  equation  for  p: 


p  •  8in(0„)  -  cos(0en,)  •  co8(0ga)  + 

+  C08(0)  •  cos(0^  •  CO8(0^^  sin(0)  •  sin(0^  -  cos(0gg^  •  cos(0Hn) 


1  -cos(0) 
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which  leads  to: 

sin(6) 

^  1  -  cos(6) 


(2.12) 


OmclusicMi 

The  vector  notations  of  the  measuring  directions  of  the  dual  lidar  system  and  the  axis  of  rotation 
of  the  auxiliary  lidar,  on  tqp  of  tiie  main  lidar.  have  been  derived.  The  azimuth  and  elevation 
angles  of  the  auxiliary  lidar  have  been  expressed  in  tiie  azimuth  and  elevation  angles  of  the  main 
lidar  and  the  angle  between  the  axes  of  the  two  systems. 
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3  WIND  VECTOR  EXPRESSED  IN  THE  SPATIAL  POSITIONS  AND  THE 

TIME  OF  FLIGHT  OF  IDEN  TIFIED  STRUCTURES 

The  principle  of  sensing  the  wind  vector  remotely  widi  an  incohemnt  hdar  is  t>ased  cm  locating  of 
moving  patterns  in  versos  range  data  s^.  Widi  die  dual  lidar,  two  of  these  sets  can  be 
recorded  simultaneously.  If  corresponding  signatures  can  be  identified  in  both  data  sets,  the  wind 
vector  <*«"  be  calculated  from  the  locations,  the  actual  time  of  the  measurements  and  the  gecanetry 
of  the  system.  The  objective  of  this  section  is  to  derive  equations  for  the  wind  vector  based  on 
those  parameters. 

It  is  assumed  that  the  structures  drift  linear  in  tiine  and  space  (during  die  measurement  period) 
along  the  wind  vector  y/,  as  illustrated  in  Hgure  3.1. 


R^neS.l:  Thiee-dimeasioiial  sketch  of  an  ttmonheiic  structure  drifting  widi  die  wind.  IhdicMed  are 

die  reference  positioo  £o(t^.  die  wind  vector  ^  and  die  position  of  die  structure  £(t)  at 
dmet 
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The  position  £  of  such  a  structure  can  be  described  in  the  time-qwce  domain  according  to: 


E(t)-^(to)  +  (t-W3K 


(3.1a) 


Writing  £(t)  ocplidtly  in  die  diree  orthononnal  directions  gives: 


E(t)  = 

X(t) 

Y(t) 

s 

Z(t) 

Po«(*o)  +  (t-^*Wco8(ogco8(oJ 

P<»y(0  +  (*-0  ‘  Wcos(oQ-8m(a,5 
+  (t-to)  •  W-sin(og 


(3.1b) 


in  which 

SB  elevation  angle  of  the  measurement 

SB  azimuth  angle  of  the  measurement 
£o  SB  arbitrary  reference  posidtm 
W  =  wind  vector 
W  »  wind  speed  (modulus  of  W) 

The  princ^le  of  the  measurement  and  the  mediod  of  die  data  analysis  has  been  visualized  in 
Figure  3.2.  The  two  lidars  measure  simuhaneously  die  atmosidieric  structures  (which  drift  with 
the  wind)  in  two  different  directions.  Each  acquired  waveform,  both  from  dw  main  lidar  and  the 
auxiliary  lidar,  represents  a  range  dependent  atmospheric  reflectioa.  The  waveforms  are  plotted 
vertically  in  false  colour  (or  in  gray  tones)  in  die  coneqxmding  figure.  Consecutive  measure¬ 
ments  are  plotted  along  the  x-axis.  If  the  rqietition  rate  of  the  lidar  is  high  enough,  the 
displacement  of  the  structures  is  uncovered  from  die  range  versos  time  figures. 
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Figure  3.2:  The  meararement  princiide  of  die  dual  lidar  and  the  method  of  prereatiiig  the  as  a 
functton  of  lange  (in  false  ctdoors  or  in  gray  tones). 


If  die  ouun  lidar  detects  a  structure  at  moment  tg,  at  range  Rg,  dim  dm  position  of  diat  Mmctme  £ 
widt  leqiect  to  an  arbitraiy  lefereace  vector  £^(l^  is: 

Em(*iii)  -  Rm  ■  mCW  •  Eo(*o)  +  (*m-<o)  ’  St  (3.2) 

in  vdiich  m  indicates  die  unit  vector  in  die  direction  of  die  measurement 

The  auxiliary  lidar  detects  diis  stmctiue  at  moment  tg  and  at  range  Rg.  This  can  be  described  as: 

E,(g  -  R.l(g  -  Eo(io)+(ta*y-3K 


(3.3) 
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in  which  a  is  the  unit  vector  in  die  measuring  dtrectum  of  die  auxiliary  lidar  at  mometit  t,. 


If  it  is  assumed  that  the  patterns  drift  linear  in  die  qiace-time  domain  (e.g.  no  rotatkm)  and  diat 
the  distortion  due  to  die  different  cross  sections  of  die  laser  beams  is  small,  die  wind  vector  can  be 
found  from: 


(t»*tin) 


(3.4) 


Substitution  of  (3.2)  and  (3.3),  while  eliminating  die  reference  position,  gives: 


= 


Ra  •<(<») 

(ta"  iffl) 


(3.5) 


Substitution  of  (2.1)  and  (2.5)  in  (3.5)  leads  to  the  solution  of  the  wind  cotiqioaents  in  the  ditee 
directions  of  the  orthononnal  system: 

Wj,  *  (jkI  •  co8(a^  •  co8(a,)  = 


R,  •  [C08(e)  •  CO8(0qn^  •  cos(0Ma4a)  »^(Q)  •  S”(0-^  «.)1 

(*a”^ 


(3.6a) 


Rju  COS(0^m,tni)  ®®®(dani,tin^ 

(ta-t«) 


Wy  »  ls(l  '  008(0,)  ■  “o(<*a)  * 

^  R,  •  [008(8)  •  co8(0qa,a)  •  sin(0,n,ta)  '  ”°(Q) '  co»(g»./.)l 

(ta-*m) 

_  Rm  •  “>»(4taMm)  “  ««°(4-m  *«) 

(ta*tin) 


(3.6b) 


w,  »  la:l  ^(og- 

^  R«  •  oos(9)  •  oos(0a^J  •  sin(0„g  -  R,,  •  8in(0«, ,») 

(«.-«») 


(3.6c) 
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where: 

|;gj|  «  W  B  modulus  of  Sf 

0  B  angte  between  dre  mein  lidtr  and  tibe  auxiliary  lidar 

^coMa  elevatknofdiemaiiilidaratdiemaiDeottit 
0^  ,.  B  azimuth  of  the  main  Udar  at  die  moment  tg, 

The  wind  qieed  Is^l  can  be  calculatcid  from: 

1^1  =V(Wjj2  +  W,2  +  W,2)  (3.7) 

The  horizontal  wind  direction  can  be  calculated  from  and  Wy  according  to: 
a,  «  2ii  -  arctan(Wy/W,)  (3.8) 

The  tenn  2ii  is  required  because  the  wind  direction  is  defined  clockwise  with  reqiect  to  die  Noidi 
vriiile  the  angles  in  the  Cartesian  coordinate  system  are  measured  counter  clockwise;  it  is  assumed 
diat  die  x-axis  is  parallel  widi  die  N<n1h. 

The  vertical  direction  can  be  calculated  by: 

aeBarc8in(W^lffil)  (3.9) 

Conclusion 

Equations  have  been  derived  to  express  die  wind  vector  in  measurable  quantities  of  identified 
structures  and  the  geometry  of  die  dual  lidar. 
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4  CAUSATION  OF  DWFTVBLCXrrrY.CHARACreRISTlC  LIFE  TIME 

AND  CHARACIERISIIC  SIZE  OF  THE  STRUCTURES  USING 
CORRELATION  TECHNIQUES 

4.1  Introdaction 

Ri  die  previous  sectioo.  a  sirai^btforwaid  m^hod  has  be»  presented  for  calculating  the  wind 
vector  in  a  plane  using  the  dual  lidar.  Apart  from  die  wind  vector  in  die  three  dimensions,  it  is 
also  possible  to  d^ennine  die  characteristic  pnqierties  of  die  structures  using  a  more  extended 
analysis.  This  is  discussed  in  diis  section. 

The  princ^le  of  inverting  die  wind  vector  from  die  drift  of  atiiioq;dieric  structures  using  an 
inodierent  lidar,  is  based  on  the  work  of  Briggs  (1950).  At  diat  dme,  die  drift  of  ionospheric 
irregularities  was  detennined  from  the  tenqxaral  behaviour  of  reflected  radio  waves,  observed  at 
three  different  locations  (three-point  measarement).  This  princqile  is  also  applied  here  to 
determine  die  wind  vector  from  lidar  measurements,  but  Iceqimg  in  mind  that  the  lidar  will 
provide  information  as  a  function  of  range  (line-sensor). 

The  subject  oi  determining  die  wind  qieed  and  die  characteristic  properties  of  the  structures  with 
lidar  in  a  one-dimensional  space  will  be  discussed  first  This  is  followed  by  a  descrqidon  of  a 
method  to  invert  the  parameters  with  two  point-sensors  in  a  one-dimoisional  ^lace  which  is 
fundamental  for  calculating  of  die  wind  vector  in  a  two-  and  in  a  three-dimensional  space. 

4.2  One-dimensional  situatioo  widi  a  line-sensor 

With  a  hdar,  a  profile  can  be  measured  along  a  inw.  within  a  fraction  of  a  second,  This  techniqoe 
has  been  dretched  in  Figure  4.1. 


TNOrapart 


Pige 

19 


— — 
■V.C' 

X 

- > 


Figiire4.1:  PtincqsleafmeasivmgpMttenu  witfalklarma(xw-dimensknalq;»ce. 

A  Mfimlatiurf  liHar  signal  is  shown  in  Hgure  4 J!a  the  attenuation  of  Ae  laser  beam  over  padi 

of  intereat  will  not  be  considered  here.  After  range-coaection  and  subtraction  of  dte  running  mean 
(called  'detrenting',  after  Stull,  1989),  onfy  the  profile  is  left  Tliis  is  shown  in  Hgure  4.2b.  First, 
Ae  characteristic  sine  and  die  diaracteristic  life  time  of  die  structures  are  derived  from  diese  data. 


Pignie4J:  Hctire  eaampie  cf  a  the  atmciapliwk  sttncture  m  measured  wkfa  lidar  (a).  The  same  signal 

afl«  detnniing  having  die  sinctuee  on^  (b). 
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Accoidiiig  to  Briggs  (19S0)  the  characteristic  size  of  die  strnctuies  is  detennined  at  die  50% 
hei^t  points  of  die  qiatial  autocooelatiaii  finctioo.  This  is  shown  in  Hgine  4.3. 


Rgure4.3:  Hcdve  esani|rie  of  the  sp«ial  autocondatioo  fiinctioa  of  a  lidar  letum.  The  characteristic 
size  of  the  stfoctnres.  X- is  determined  at  die  50%  pomts  of  the  maximum  value  (Rill  Width 
HalfMaxiinmn). 

If  the  structures  vary  in  dme  only  (no  drift)  than  the  value  of  die  cross-correlation 

function  of  two  successive  records  (same  padis  but  different  times)  decreases  with  increasiag  time 
between  die  measurements.  According  to  Briggs  (1S150),  the  diaracteristic  life  time  is  equal  to  die 
time  delay  between  two  successive  records  of  which  die  (tenqioral)  crossKXsrelation  reduces  to 
50%.  Hgure  4.4  shows  a  simulation  of  die  variation  of  die  cross-correlatitm  as  a  ftmction  of  the 
time  delay  b^ween  two  records  and  the  princ^le  to  detennine  die  characteristic  life-time,  t^ 
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Figine  4.4:  Sinrmtated  variation  of  tbe  temporal  crosa-^onelatioa  venus  tbe  tune  deli^  between  two 
lidar  records.  Ibe  characteristu;  life  time,  t^  is  defined  as  diat  time  d^y  where  the 
maximum  value  of  the  temporal  ooss-coirelatioa  has  been  decreased  to  50%  of  the 
maximum  value. 

Thus  far,  only  die  spatial  a^tocoIte^tig^l  function  of  one  record  is  considered  and  die  temporal 
cross-correlation  function  between  successive  records.  But,  if  die  tune  series  at  range  'x'  is  cross- 
correlated  with  tbe  time  series  (same  time  interval)  at  x-fAx,  dian  a  point  of  die  cross- 

correlation  function  is  obtained.  By  varying  Ax  over  die  whole  sptct,  die  qmtial  cross-correlation 
function  is  obtained.  Because  die  structures  will  also  change  in  tune,  die  maximum  value  of  the 
spatial  cross-correlatum  function  decreases  as  die  time  increases.  Spatial  cross-correlation 
functions,  for  diffoent  time  intervals  but  continuous  in  qmce,  are  shown  in  Hgnre  4.5a.  A  two- 
dimensional  correlogiam  can  now  be  constructed  from  a  large  set  of  lidar  records  as  shown  in 
Rgnie  4  Jb  (being  die  top  view  of  figure  4.5a  where  die  ptdnts  of  equal-correlation  are  cmmected 
by  dashed  lines).  Hgure  4  Jc  show  die  same  correlogram  if  space  and  time  can  be  varied 
continuously. 
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In  dte  next  stq>,  we  diet  die  atnictnres  in  our  one-daneiuional  qiace  will  also  drift  widi  a 

constant  velocity.  This  means  that  the  temporal  variations  in  the  structures,  as  measuied  with  (be 
lidar,  are  caused  by  two  independent  e£Eects:  drift  and  turbulence.  Due  to  die  drift,  die  (peak  of 
the)  ^Mtial  cross-conelatkm  function  of  two  lecosds  is  shifted  scmoevdiat  in  ^Mce,  as  shown  in 
Hgure  4.6a  (time  axis  petpendicuiar  to  die  figure).  The  shift  is  equal  to  the  product  of  &e  drift 
veloci^  and  the  time  del^  between  die  measurements.  Due  to  turbulence,  die  maximum  value  of 
the  CTOSS-coiTelation  function  has  been  decreased.  If  the  spatial  cross-conelation  fiinctions  are 
plotted  in  vertical  planes,  sqiarated  by  time  intervals  of  At,  again  a  quasi  doee-dimensional 
coirelogram  arises.  The  top  view  or  die  projection  of  that  correlogtam,  in  which  the  points  of 
equal  cocrelatkm  are  connected  by  dashed  lines,  is  shown  in  Figure  4.6b.  Hgure  4.6c.  shows  the 
projection  of  coirelogram  if  the  time  can  be  changed  continuously. 

Due  to  die  drift  of  die  structures,  the  axes  the  of  the  iso-cross-coirelatitm  ellqises  are  rotated 
someudiat  as  can  be  seen  in  Figure  4.6b  and  4.6c.  Nevertheless,  it  is  still  possible  to  determine  die 
characteristic  life  time,  tc.  die  characteristic  size,  x^,  and  die  velodQr,  V,  of  the  drift  from  this 
figure  as  will  be  shown  in  the  next  secdixi. 

If  there  was  no  drift,  die  characteristic  life  time  and  die  characteristic  size  were  determined  frxim 
the  inlercqrts  of  die  elli^,  with  a  cotrelation  value  of  0.5,  widi  respectively  die  vertical  and  the 
horizontal  axes.  But  in  die  presence  of  drift,  the  initial  inteicqit  of  die  ellipse  widi  the  vertical  is 
shifted  in  die  horizontal  direction  by  an  amount  equal  to  die  product  of  die  drift  velocity  and  die 
time  delay.  This  means  diat  in  the  presence  of  drift  the  characteristic  life  time,  tc,  is  d^ennined  by 
die  horizontal  tangent  to  the  eUqise  with  correlation  value  of  0.5.  Because  there  is  no  shift  in 
horizontal  direction  <m  die  x-axis  (Ml),  die  characteristic  size,  is  given  by  the  intercqit  of  diis 
elli^  with  die  hmizmital  axis. 

The  drift  velociqr,  V,  can  be  determiiied  fimn  the  ratio  of  die  characteristic  size  and  die  diaracter- 
istk  life  time  but  it  is  also  possible  to  use  dw  horizontal  intercqit  ai^  odier  ellipse  and  its 
vertical  tangent 

rmirinMnw  fhe  two-dimensional  time-qpaoe  correlogiam,  as  derived  from  a  set  of  consecutive 
lidar  measurements  in  a  one-dimensional  qpace,  provides  both  die  duttacteristic  properties  of  die 
structures  as  wdl  as  their  velocity.  The  tedudqne  diarnff^  boe  can  be  applied  to  determine  die 
radial  wind  in  the  horizontal  {dane. 


TNOr«port 


Rgiiie4.6: 
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on^  axis  is propoilioiial  wWi the tfaie 

intoval hnrw!!nrtr?*  •  of  die  two-dimenido^  condofiam  if  die  tfaoe 

laiervaipetween  the  wcoesshreiecoids  can  be  inade  infinitely  small  (cjr^ 
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4.3  One-dimensioiial  sitiuttioa  widi  two  point-sensors 

With  two  point-sensors,  it  is  also  possible  to  derive  the  characteristic  life  time,  die  diaractetistic 
size  and  the  velocity  of  die  stroctme.  This  situation  occurs  e.g.  when  the  dual  hdar  is  used  for 
cross-wind  measurement  as  a  function  of  altitude,  as  shown  in  Figure  4.7.  If  die  wind  direction  is 
known,  than  the  path  between  the  aauon  (sait^ltng  points)  is  defined. 


DUAL  LIDAR 


Figure  4.7:  Schematic  set  up  for  measuring  the  cross-wind  speed  as  a  functioa  of  height  widi  a  dual  Udar 

system. 

At  die  moment,  die  objective  is  to  determine  die  wind  qieed,  given  the  wind  direction  and  die 
distance  between  die  point-sensors,  xq.  Random  stmctmes  move  along  die  one-dimensional  padi 
and  are  recorded  as  time  series  from  two  different  locatkmt.  For  die  determination  of  die  wind 
qieed,  die  time  series  from  diose  two  locations  are  croas-coaelated.  The  quotient  of  die  distance 
between  die  sensing  points  and  die  time  shift  where  the  maximum  cross-corrdatkm  occurs, 
provides  the  drift  veloci^. 

The  determinatiQn  of  die  dtaracterisdc  size  and  die  characteristic  life  time  is  somewhat  different 
From  die  point  of  view  of  die  two-dimensional  cross-coirelogram,  as  diown  in  Figure  4.6,  onty 
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two  twiiporal  autoconelations  are  available  (which  should  be  ulmtical  if  the  structures  can  be 
described  with  an  ergodic  process)  and  one  tenqxaal  cross-conrelaticMi  function.  This  has  been 
sketched  in  Figure  4.8. 


Figure  4.8:  Autoconelatkm  and  cross-coneliuion  functions  of  two  time  series  measured  at  two  different 
locations  (a)  and  their  positioiu  in  the  two-dimensional  conelogram  (b). 

The  inversion  of  the  drift  velocity  and  the  characteristic  properties  of  the  structures  is  straight 
forward  but  requires  extensive  maOtematical  evaluatioiL  Therefore  only  the  necessary  basics  of 
the  solutions  will  be  discussed  here.  The  full  derivation  of  the  two  different  methods  have  been 
elaborated  in  the  Appendices  A,  B  and  C. 

The  analysis  is  based  on  die  assunqition  diat  the  iso-conrelation  lines  can  be  described  by  ellqises 
(see  Figure  A5  and  4.6)  as  discussed  by  Briggs  (1950)  and  Zuev  et  al.  (1977).  As  shown  in  Hgure 
4.8b,  die  autocorrelation  fimctirms  of  the  tune  series  measured  with  die  two  point-sensors  are 
plotted  in  the  vertical  plane  at  x^O  (die  results  should  be  identical  for  the  two  time  series).  The 
cross-conrelation  function  is  plotted  in  die  vertical  plane  at  x^xg  whidi  is  equivalent  with  die 
distance  between  die  sensors.  The  objective  is  to  derive  both  die  characteristic  size  and  the 
characteristic  life  tune  of  the  structures  as  weU  as  the  drift  velocity.  This  is  done  by  calculating 
die  length  of  the  axes  of  the  ellipse  and  die  angle  of  rotation. 
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Mediod  1: 

A  Active  iso-ooss-cofielatioa  ellqise,  constnicted  fiom  the  antoconelatioo  function  and  die  cross- 
conelatioo  function  as  shown  in  Figure  4.8b,  provides  the  cooedinates  (xq^q)  and  (O^^).  A  lengdiy 
algebra,  whidi  has  been  elaborated  in  .^ipendix  B,  provides  die  following  results  eaqiressed  in 
direct  measurable  quantities: 


The  characteristic  life  time: 

tci  -  V(to2+ti2) 

The  characteristic  size: 

_ 

***  ■  >/(to2+ti2) 

The  wind  qieed: 

V  - 

‘  (to^ti2) 


(4.1) 


(4.2) 


(4.3) 


Method  2: 

If  it  is  difiicult  to  detennine  die  maximum  value  of  die  cioss-conelation  function  than  it  is  also 
possible  to  use  die  value  of  diis  function  at  the  time  shift  tj  where  this  function  becomes 
again  equal  to  die  function  value  at  and  the  time  shift  t2  where  autoconrelatioo  functitm  value 
(0.3)  is  equal  to  die  cross.<oiielati(xi  value  fimetion  at  M).  With  these  diree  coordicates  one  can 
also  calculate  die  characteristic  size,  die  characteristic  life  time  and  the  drift  velocity.  The  algebra 
has  been  elaborated  in  Appendix  B  and  die  results  ate  given  in  (4.4)  to  (4.6) 


The  characteristic  life  time: 


*c2 


(4.4) 


The  characteristic  size: 
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The  velocity  of  (be  drift: 


tl  -xl 


(4.6) 


CoftchiMon 

Id  dw  ooe-dimenaioaal  atnarioo,  die  wind  speed  and  the  characteristic  ptopecties  of  the  eddies 
can  be  derived  fiom  die  dme  series  of  the  data  lecosded  widi  two  point-aeiisars  some  distance 
tpait  The  objective  is  reached  by  analysing  die  auto-  and  cross-conelation  functions  obtained 
ftom  leconled  data  sets.  The  only  requireinent  is  (hat  the  time  series  are  long  enough  to  describe 
die  patterns  widi  sufBcient  statistical  reliability. 


4.4  Two-dimensional  situation  widi  two  line-senson 


The  horizontal  wind  speed  and  wind  direction  can  be  inverted  frun  two  sets  of  lidar 
measurements  reconled  simultaneously  in  two  different  direcdcms  as  shown  in  Figure  4.9. 


ngiiie4.9:  Set  im  for  measuring  the  horizomal  wind  vector.  The  dme  series  at  point  R- is  cross 
comhaed  with  an  die  time  series  of  the  jpofatts  along  padi  a.  The  time  sei&s  with  the  laq^ 
ooss-oonelation  vahe  determines  poshion  whOe  die  time  riiift  rrf  that  mazhiinm  value 
determines  die  transit  time. 
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The  calculation  stam  with  the  tisne  aeties  of  the  atmoqpheric  hackacattca'  at  nnge  This  time 
series  is  cross^onelated  with  all  availaMe  time  series  as  measured  with  die  second  lidar.  Rxxn 
die  set  of  cross-coadatioD  fhnctioiis.  die  one  udiidi  provides  die  largest  cross-canelatimi  value 
indicates  the  positian  The  time  delay,  of  the  maximum  value  in  the  cooeqwnding  cross- 
correlation  function,  is  equtvaknt  widi  die  twie  required  for  the  structures  to  move  from  posidra 
Rg,  to  position  R,.  As  a  result,  die  wind  vector  can  be  calculated.  This  princ^k  is  equivalent  widi 
the  mediod  described  in  previons  secdcm  and  has  been  used  to  derive  die  horizontal  wind  vector 
as  described  in  section  S.2. 

4.S  Three-dimensional  situations  widi  duee  line-sensors 

The  wind  vector  in  a  horizontal  plane  can  be  derived  from  horizontal  lidar  measurements 
performed  in  two  difierent  directions  as  shown  in  section  4.4.  By  expanding  diis  idea  to  duee 
lidars,  the  wind  vector  in  duee  dimensions  can  also  be  found.  (It  is  assumed,  diat  widiin  die 
dimensions  of  die  eiqieriment,  die  aerosol  eddies  drift  only  in  a  linear  direction  without  rotation.). 
The  practical  realization  of  this  princ^k  has  been  sketched  in  Hgure  4.10. 


Rgure  4.10:  Ftinc^  of  measuring  the  duee-diiiiensional  wind  vector.  The  sani|riing  poktt  A,  B  and  C 
are  not  necessary  in  a  horizoaial  pbaie.  Locttkn  A  detenninet  die  reqjiiired  attitude.  The 
t  tune  scries  measured  at  A  is  cross  correlated  with  all  available  time  series  along  die  ptfis  2 

and  3.  Those  thne  series  rriikh  provide  die  maximum  cross-oonelation  values  mermiiie  die 
locadoiu  B  and  C. 


bi  tfiis  caie,  ttae  time  series  of  die  siinoi^iheric  sttiictiBe,  mrs sored  widi  die  first  lidv  (A)  st  s 
selected  locstion  Qiei^)  is  cross  ccnelsted  wtdi  Ike  time  series  of  Ike  second  lidsir  (2)  over  sU 
svailaUe  locstioos  (heig^).  l¥om  due  set  of  cross-oonelsrion  fimcrimps  Ike  one  diet  provides  the 
largest  ctossKxxielalioii  vahie  is  dwsen  to  belong  to  Ike  most  likety  second  location  (B).  The 
time  riiift.  vriiich  coneqionds  to  die  mazimum  valiie  in  die  cross-coodation  fimctioo,  is 
equivaknt  vridi  the  time  required  for  Ike  structures  to  drift  ftom  position  (A)  to  positum  (B).  In 
die  same  wi^,  Ike  best  condatkm  of  die  time  series  ftom  position  (A)  is  selected  from  die  time 
series  of  lidar  3.  This  provides  die  diird  location  (Q  and  die  second  transit  time.  The  dnee 
locations  A,  B,  and  C  define  the  plane  of  die  wind  vector  »««nmm£  die  structures  drift  widi 
the  wind  and  diat  die  size  of  die  structures  and  die  selected  locations  do  not  influence  the  results. 
The  triangle  ABC  has  been  redrawn  in  Rgi  4.11. 


ngnn4.11:  Geometry  of  the  three  snnpling  locations  as  defined  in  Figure  4.10  from  whkk  the  wind 
vector  is  calculated. 


Calculation  of  die  wind  vector 

For  large  structures,  die  transit  time,  At^,  to  drift  from  A  (thus  from  ^  to  B  is  determined  by  die 
ratio  of  die  path  lengdi  (?B  and  die  wind  qieed.  hi  die  same  wi^,  die  time  required  to  drift  from 
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A  (from  0)  to  C,  Atj,  equals  die  ratio  of  QC  and  die  wind  qpeed.  To  derive  die  wind  vecttH',  first 
die  angles  a  and  ^  are  calculated  i^iich  define  the  udiole  geom^  of  die  figure.  Then  the  wind 
direction  is  givm  by  die  unit  vector  along  Q3  (or  along  QQ  and  die  wmd  can  be 
calculated  from  the  ratio  of  Q'B  and  At2  (or  the  nuio  of  and  Atj). 


Due  to  the  measuring  principle,  the  wind  vector  £  is  in  the  shaded  plane  ^  which  will  be 
described  by  the  vectors  It  and  £  being  re^ectively  die  unit  vectors  akmg  AC  and  AB.  If  i  is  a 
unit  vector  akmg  Q-Q*  than  the  scalar  products  of  I  widi  b  and  with  £  provide: 

I .  It  »  cos(180-<x)  ■  -cos(a)  (4.7) 

I .  £  »  cos(P)  (4.8) 

Because: 

Q'B  »  l3Kl  Atj  «  AB  sin(P)  (4.9) 


QC  «  Iwl  Atj  *  AC  sin(a) 


(4.10) 


one  can  eliminate  |£|  and  find  the  relation  between  a  and  p  eiqpressed  in  measurable  quantities. 
This  results  in: 


AB  ■  sin(P)  AC  -  8in(a) 
Atj  At} 

Because: 

o  +  p  +  y«  180 


(4.11) 


(4.12) 


in  uduch  yis  givenby  die  getmwtty  of  die  nqieriment  and  die  locations  of  A.  B  and  C,  a  second 
equation  is  available  to  solve  for  both  angles.  Substitution  of  equation  (4.12)  in  (4.11)  and 
eliminating  a  leads  to: 

At}  AC  8in(y)  •  cos(P)  +  sin(p)  •  cos(y) 

At2  AB  sin(p) 


(4.13) 
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This  provides: 
tan(P) 


B»(10 


Mx  AB 


(4.14) 


coidi 


Atj  AC 

KKtnfiMrinn  of  ^  from  equstioiis  (4.11)  and  (4.12)  leads  in  die  same  way  to  a  aohition  of  a. 
tmdi 


tan(a) 


Ah  AC 

12. - C08(Y) 

At3  AB 


(4.15) 


The  wind  speed  can  be  calcnlated  ehher  from  equation  (4.9)  or  from  equation  (4.10)  using  eidm 
equation  (4.14)  or  (4.1S)  reflectively. 


Expression  &g  die  wind  ditecticm 

An  auxiliaiy  vector  b,  vrtiich  is  perpendicttlar  to  plane  can  be  described  by  die  unit  vectors  b 
and  £  According  to: 

b  «  bx£  (4.16) 

Note  that  die  cross  product  of  die  unit  vectors  i  and  b  provides  die  unit  vector  in  die  wind 
diiectkm  (s><ixb)  vriuch  can  also  be  described  by  die  triple  product: 

ixb-cx(bx£)  (4.17) 

The  left  hand  ride  of  (4.17)  is  equal  to  die  unit  vector  in  die  wind  £.  Elaboration  of  die  ri^t  hand 
side  ci  (4.17)  results  in: 

S  -  b([£)-£(X.b)  (4.18) 

Subsdtidion  of  equations  (4.7)  and  (4.8)  provides  die  unit  vector  in  die  wind  diiectkin: 

X  -  b  ‘Cos(P)-t-£  oos(a)  (4.19) 

GoBCbiiyflB 

A  mediod  has  been  described  to  determine  the  duee-dimensional  wind  vector.  This  mediod  is  one 
stqi  fufdier  dian  diose  publidied  so  fiv,  which  derive  only  die  horizontal  wind  f.«wipnnt»n*  as  a 
ftmcdon  of  ahftode.  The  new  mediod  ndghi  be  interesting  for  mfiping  convection.  The 
set  of  equations  derived  in  diis  section  and  in  (he  prevkms  section  win  be  used  in  die  next  sectkn. 
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S  examples  ex' UDAR  RESULTS  AND  THE  DERIVED  WIND  VECTOR 

5.1  Tntrodnetioo 

RTpiaimenta  OH  the  detenoiiuttioo  of  the  wind  vector  wifh  tibe  Mie  hdir  started  in  the  qsing  of 
1991.  After  a  period  <tfdevek^«ig  software,  oq^ioviag  dte  laser  and  improving  die  receiver,  the 
routme  measaremeats  were  started  in  Angnst  1991.  IhMe  measnreaaeats  were  performed  in  a 
fixed  hnriaiiitoi  and  fixed  vertical  direction  to  detmnine  die  temporal  bdutviour  of  die 
atmosphere.  The  lidar  was  used  also  in  die  scanning  mode,  in  a  horizontal  and  a  vertical  plane,  to 
m«p  the  atmospheric  strectore.  In  a  later  stage,  die  triangulatkm  measurements  were 

performed  as  well  to  detennme  die  wind  vector  as  a  function  of  height 

Some  of  die  data  sets  obtained  with  die  dual  lidar  are  shown  in  this  cluqiter.  Section  5.2  describes 
a  subjective  wnpHMutinn  (rf'  die  amount  of  structure  generally  found  in  die  data  sets  (the  structure  is 
necessary  for  deriving  die  wind  vector).  In  section  S.3,  some  results  of  the  horizontal 
measurements  ate  presented  and  in  section  5.4,  some  results  of  die  vertical  measurements.  In 
section  5.5,  some  results  are  jnesented,  obtained  from  triangulation  measurements  to  a  maximum 
altitude  of  1000  m. 

5.2  Statisticaloverviewof  the  results  in  the  period  Sqitember  1991 -January  1992 

In  die  period  September  1991  -  January  1992, 82  successful  eiqieriments  were  carried  out  widi  die 
dual  lidar.  These  experiments  consist  of  a  set  of  horizontal  and  a  set  of  vertical  measurements. 
The  structure  in  die  data  sets  has  been  visualized  by  subtracting  die  average  vahre  from  each 
sample  and  by  du|di^ing  in  fidse  colonr  only  die  variations  within  duee  times  the  standard 
deviation.  The  structures  obtained  were  subjectively  cL.ssified  ftom  0  (no  structure)  to  5  (much 
structure).  Table  5.1  shows  the  results  of  dus  arudysis. 
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TableS.l:  Nanoiier  of  dtta  nis  aad  dieir  pcrceataget  in  nfajecttve  ^uet,  rapmaitiag  die  ainoint  of 

imictaie.  riowet  4  and  5  indi^  due  die  due  ect bepoueeiid  candidue  tor  inveniao  of 
the  wild  vector.  Ibe  lidar  awaaonneBea  irae  performed  b  dw  period  Scptendier  1991  to 
jamaiy  1992. 


aaas 

W.W  — • - 

rnuTWigr 

Horz.  Vert 

Percentage 

Horz.  Vert. 

0 

49 

12 

60 

15 

1 

16 

30 

20 

38 

2 

6 

15 

7 

19 

3 

1 

7 

1 

9 

4 

6 

7 

7 

9 

5 

4 

8 

5 

10 

Total 

82 

79 

100 

100 

The  reaiilta  of  TaUe  5.1  indicate  diat  in  only  12%  (dass  4  and  S)  of  the  cases  tbeie  is  suffickot 
structure  to  detennine  die  hodzoatal  wmd  vector.  In  die  vertical  direcdcMi,  19%  of  the  data  sets 
diow  sufiBcient  structure. 

The  large  number  of  occasknis  widi  insufficient  atmoqiheric  structure  might  be  caused  by  die 
perfect  mixing  process  of  die  aerosol  eddies  at  low  altitude,  ff  inversion  layers  are  present,  modi 
structure  is  observed  at  die  top  of  die  boundary  (entrainment  layer),  but  diis  interval  is  too  limited 
to  detennine  a  wind  vector  within  die  whole  mixed  liyer.  Above  die  mixed  layer,  however,  where 
less  aerosols  are  present,  no  structures  have  beoi  detected. 

To  d^etmine  die  wind  vector  as  a  fbnctkn  of  altitude,  die  lidar  was  pointed  in  three  different 
(elevated)  dhectioiis  as  described  in  section  4.5.  In  the  period  November  1991  to  January  1992, 
37  triangulation  measurements  were  petfonoed.  The  dedaion  to  perform  such  a  complicated 
measurement  was  based  on  results  of  die  vertical  measuremenL  The  data  sets  were  classified  as 
shown  in  Table  5.2. 
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Table  5.2:  PialiawtkiB  of  the  daiaea  ot  the  rewltt  of  ifae  triangwlalioa  aeamtcmeBtt.  See  alio  text 

wMi  table  S.l. 


C3sss 

Number 

Fetcentage 

0 

4 

10 

1 

10 

27 

2 

10 

27 

3 

8 

22 

4 

4 

11 

5 

1 

3 

Total  37 

100 

Only  die  clasaes  4  and  5  provide  sufBcieot  infonnatum  to  detexmme  the  wind  vector  from  die  «<«*« 
sets.  This  means  only  in  14%  of  the  cases,  whereas  this  ^pe  of  measurement  was  only  performed 
if  the  vertical  measurement  provided  enoo^  information.  The  reason  for  this  low  percentage 
might  caused  by  the  relative  long  time  interval  for  diis  type  of  measurement 

5.3  Some  results  frwn  horizontal  measurements 

Some  results  obtained  from  horizontal  measurements  are  presented  in  this  section.  For  clarity, 
first  some  results  for  cases  widi  no  atmoqiheric  structure  are  presented.  Next,  some  sets  are 
shown  widi  sufficient  structure.  At  die  end  of  diis  section,  an  overview  is  given  of  the  horizcmtal 
wind  vector  determined  from  anumber  of  data  sets. 

Data  sets  with  no  stroctore 

Two  mtamples  of  data  sets  udiidi  show  no  or  very  little  structure,  are  presented  in  Rgmes  5.1  and 
5.2.  The  actual  meteorological  sitnatkm  is  printed  in  the  figures. 

The  vertical  lines  in  the  figures  ndg^t  be  caused  by  variation  of  die  reference  level  or  die  laser 
energy.  In  diis  stage  of  die  experiment  it  is  not  possible  to  compensate  for  diese  ^Gects.  The 
horizontal  lines  are  on^preaent  in  die  fignuM  and  do  nnrinfliMnr<>iV!  film] 

The  lidar  data  aet  shown  in  ngure  5.2  has  acme  stroctnre  in  the  first  half  of  die  period  but  widi 
limited  tignal-ioHiidae  ratio. 
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FOe:  911123-240Q29-1630 


Mfiloa: 


Tempetaluie 

-  M 

±  Ofi 

•c 

-  88,7 

±  3A 

% 

Wet  bulb 

- 

±  Ofi 

•c 

Sdntaiatkn 

-  102 

±  Ofi 

- 

ReL  HumidiQr 

-93^ 

±  0.4 

% 

Waidipeed 

-  2.65 

±  0,7 

m/f 

Fnaanre 

-  1009,7 

±  0,25 

hPa 

Winddfaecdoa 

-  142,2 

±  3,9 

Degr. 

\^aibiliqr 

±  OfiS 

km 

FtgoieS.l:  Range  Ycwntime  fiyet  of  die  Mmnmhrric  iciWwring  ai  meMUWd  wilfa  two  IMmi.  I^yer 

Rgnm*  Raw  data  aet  aa  recoaded  win  die  Rdan.  Lower  figona:  Same  data  aet  aa  qm 
figmea  after  of  die  rowa.  HorixoHtal  aoda:  0-272  a;  vertkal  azia:  0-2352  m,  for 

eadiBdar. 
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Tempctamre  -  937  ±  0^  *C 

Wetbulb  -  8,00  ±  0,00  ‘C 

ReLKnnidby  -  83,71  ±  0338  % 

Ftamm  ■  102039  ±  036  liPk 

-  2,73  ±  0,04  km 


Tnumiarion  -  67,74 

Sdntflktka  -  2,46 

Windipeed  -  032 

IKMdirectioo  -  98,82 


±  338  « 

±  035  - 
±037  m/s 
±  135  Degr, 


Figiin53: 


Kaige  TCBns  tiow  fignre  of  tiiB  simoiplieiic  ■r.ttwrinj  Upper  fisuR 
noonied  Iqr  ihe  Udir.  Lower  fignie:  Sim  <hm  set  m  npper  qrai 

Horimolri  axis:  0-271  s;  vertical  axis:  0-2016  m,  ftr  eadi  lidK 


Ryv  data  set  as 
after  detraitieg. 
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Data  acts  with  stracture 

fia  some  sifwarions  die  lidar  delects  mudi  stmoqiheric  stractme.  This  may  occur  bodi  during 
(light)  precipitation  but  also  during  very  clear  atmospheric  conditions. 

The  first  exaaqile  of  an  atmoqiheie  withmneh  structure  is  shown  in  Rgure  S.3. 

Hgure  S.3  shows  the  qiatial-tenfxual  bduviour  of  die  atmospheric  aerosol  stnictures  as 
measured  widi  one  lidar.  This  example  with  so  much  structure  is  very  uncommoiL  Unfortunately, 
no  results  of  die  seoMid  Udar  are  available.  This  data  file  has  been  elaborated  in  more  detail  in 
Kunz,  1993. 

Other  exan^les  of  inhomogeneous  horizontal  atmosfdieres,  as  measured  simultaneously  with  two 
Udars.  are  shown  in  Rgures  5.4  and  5.5.  Because  the  wind  is  almost  in  die  direction  of  the  lidar, 
die  patterns  have  negative  slopes  with  reflect  to  die  horizontal  axis,  Le.  die  structures  are 
advected  toward  die  lidar  and  get  closer  betwem  each  sounding. 
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MbIbo: 

TempenlufB 

•  11,6 

±  0,02  *0 

Tnomiiasioo 

-  913 

±  1.4 

% 

Wet  bulb 

- 

•c 

Sdnlfllatiao 

-  53 

±  038 

. 

ML  iBiiniiHiy 

-  86,6 

±03  % 

Wlndqpeed 

-  531 

±  132 

Bmme 

-  994.1 

±03  bPik 

wm  QiwcnoD 

-  2053 

±  533 

Degr. 

^^■ibili^ 

-  5,9 

±0,1  km 

FigmeSJ: 

Itawr  figoie:  nw  dM  as  meuimd  by  Ae  Udv.  Ixnrar  figne:  w 
Hoiizcattd  isiK  9>128  s;  veRical  axis:  0-1392  m. 

me  data  after  detSBoting. 
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Meteo: 
Tempename 
Wet  bulb 
Rd.Ibiiiidity 
Fnaanie 


-  10,07  ±  0,10  ’C 

-  9,27  ±  0,10  ‘C 

-  95,80  ±  1,36  % 

-  1024,50  ±  0,23  hPa 

-  43,25  ±  0,004  km 


Transmiwkn  -  91^ 

Scintillation  -  2,889 

Windowed  -  2,68 

Wlnddinctioa  -  281,83 


±2^  % 

±  0,269  - 
±  0,73  m/s 
±  2,49  D^. 


Hgnre  5.4:  U 


n:  law  data  aa  meaaond  by  die  lidar.  Lower  figure:  same  data  after  delieatmg. 
ads:  0-272  a;  vertical  axis:  0-1680  m,  for  each  Udar.  See  also  figure  5.5. 
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HgineS^: 


CtndiiDed  from  figae  5 A  Upper  figure:  aw  dae  n  meaamed  by  dw  lida.  Lowa  figme: 
nme  dae  after  dettenting.  Honzonal  udK  (V272  s;  vatkel  aodK  0-16M  m,  fiv  each  leeda. 
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Afiahwii.  of  tfi*  KAmr  <!«»«  FVgVfft  Wj  S,S 

a.  Radial  wind  qieed  derived  from  die  dvectko  of  tbe  pattens 

The  slope  of  the  stmctures,  with  fcqwct  to  die  hoiizootal  axis,  provides  die  radial  wind 
speed.  The  esthnated  vahie  for  die  slope  is  about  7.06  tn/s  (1920  taP12  s)  while  die  actual 
wind  speed  was  about  2.68  ±  0.73  oo/s.  This  large  differeiice  will  be  discussed  later. 

b.  Wind  vector  determined  by  cross  correlating  time  fimctioiis 

The  dual  lidar  provides  die  possibility  of  measuring  die  transit  time  of  the  atmoqihetic 
aerosol  eddies  betwem  the  two  hdar  axes  as  discussed  in  die  previous  sections.  In  diis 
exan^le,  the  time  series  at  a  certain  range  from  die  first  lidar  is  cross  correlated  widi  die 
time  series  at  different  ranges  fiom  die  secrnid  lidar.  That  pair  of  time  soies  whidi  provides 
the  largest  value  in  die  cto8S<cotrelatk«  function,  indicates  dM  maximiiin  likelihood  of  the 
structure.  The  time  shift  of  die  peak  in  die  correlation  fimction  gives  the  transit  time.  In  diis 
exan^le,  die  time  series  of  the  auxiliary  lidar  at  750  m  is  cross  correlated  with  die  time 
series  of  the  main  lidar  between  540  and  1140  m.  The  maximum  value,  the  time  shift  and  die 
uncertainty  (90%  value)  in  die  time  diift  are  shown  in  Table  5.3. 

TableS.3:  Maximum  valuet  of  the  cross-correlatioo  function  and  the  time  duftaPani  file  24083640. 


The  results  of  die  analysis,  eiqxessed  in  actual  range  and  time,  ate  diown  in  Hgure  5.6. 
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MAX.  CROSS  CORRELATION  VALUE 
and  TRANSIT  TIME 

-  CfS-Crr  *  Shift 


RANGE  2-nd  lidar  (m) 
(RANGE  1-(t  LIOAR:  750  m) 


Rguit  5.6:  Remits  from  die  cross-condatioo  cakulsdon.  The  dots  indicste  the  maTimnm  value  of  the 
Cfoss-oonelatiaa  fuoctians  (left  scale).  Ihe  plus  signs  and  the  enor  ban  indicate  die  rime 
shift  (right  scale). 


The  time  series  fircnn  the  auxiliaiy  lidv  at  7S0  m  and  from  die  main  Udar  at  908  m  ate  shown  in 
Hguie  S.7a.  Antocoxtelation  functions  derived  from  die  main  Bdar  data,  calcnlated  at  centre  of  the 
time  series  over  different  time  windows,  are  shown  in  Kgme  S.Tb.  The  cross-correlation  function 
of  the  time  series  from  die  main  lidar  widi  die  time  series  from  the  auxiliary  lidar,  are  shown  in 
Rgure  S.7c. 


TIME  (**c) 


TIME  SHIFT  (Mc) 


TIME  SHIFT  (MC) 


i.7:  Tune  leriei  fkvuu  Figure  5^  nd  53  (u).  Autoconelrekio  ftnctioas  of  die  «Min  Udar,  at 
906  DB  mge,  in  the  ficwire  ot  (be  tine  wties  with  time  tine  windows  (b).  Ooss-ooiielatian 
taocticos  witii  die  wnOiaiy  Ikkr  n  7S0  n  noge  (c). 
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The  diaracteristic  life  time,  ^  and  die  cfavacteriatic  tiae,  die  Mractures  have  been 
painiitwH  frooi  dw  aiito-  and  oroaa-conelatiaci  ftmctions  according  to  equations  (4.1)  utd  (4J2) 
resulting  in:  23.710.5  s  nd  1845±1  m.  The  uncertainty  shown  fidlows  from  dw  variation 
(ji  die  diree  correlation  ftmctions,  however,  the  actual  uncertainty  can  only  be  estimated  after 
anatysing  all  the  time  series. 

Prom  Rgure  5.6,  die  uncertainty  in  die  position  of  die  peak  in  the  aoss^andation  function  (90% 
values)  has  been  indicated.  This  means  that  die  data  record  of  die  second  Udar  comes  from  range 
90730±17.37  m  and  has  a  delay  of  24.11±6.8  s.  Taking  into  account  die  geometry  of  the  dual 
lidar,  die  calculated  wind  speed  becomes  7.11±1  m/s  and  die  wind  direction  becomes  314.9012.2 
degrees  widi  reqiect  to  NotdL  The  actual  wind  direction  was  281.8123  degrees  and  die  wind 
speed  was  2.6810.73  m/s. 

To  find  a  reascm  for  the  discrqpancy  between  the  lidar  and  die  cap  anenunneter,  die  data  file  was 
also  analysed  at  other  locations  in  die  record.  Moreover,  odier  data  files  widi  sufficient  structure 
were  analysed.  If  possible,  also  die  wind  qieed  was  derived  firom  die  direction  of  die  patterns  in 
the  false  colour  figures.  The  results  of  diese  analyses  are  summarized  in  Table  5.4. 


Table  5.4:  Summaiyoftheresuhsof  die  wind  qieed  analyses 


dne 

Ok 

wlDdcep 

mtammrnK 

dope 

wind  derived  wtt  lidv 

index 

911024 

24063640 

248±2.49 

281J313J3 

706 

60211437 

294.191  90S 

SO 

6341  136 

292001  203 

73 

6.76il  133 

284.731  107 

100 

4.701  130 

289031  203 

123 

3041  233 

278.4911139 

130 

911023 

23100943 

0.32±1J8 

109.14±2.1S 

330 

3301  030 

116001000 

30 

3.131  04S 

118.191033 

73 

3301  048 

116381330 

73 

4301  04S 

118.191033 

73 

3071  - 

108091  - 

100 

3301  - 

102.761  • 

100 

3341  - 

113071  - 

100 

911106 

06093624 

3.4S±1J1 

2443012.7 

9.1211400 

26906133.70 

40 

73911100 

2443813800 

30 

3001302 

2340713831 

80 

91110S 

0S04S383 

3.9adbl.ll 

331.411013 

730 

6071333 

329.9911034 

40 

6071730 

690313401 

40 

3311237 

3290911034 

40 

6021308 

3310711900 

60 

7.181233 

340331831 

80 

911217 

1714S700 

Sj01±1^ 

209331533 

- 

80111433 

2763116009 

60 

911204 

04163230 

03611131 

241j0«12O1 

• 

2.741637 

269061106.0 

40 

3.741 1.18 

299301108 

40 

1.181031 

269301334 

40 

911211 

11143330 

lJ4dOj64 

msmsa 

2J0 

2.781130 

132061331 

40 

911220 

20142702 

io.io±im 

296.761136 

_ 

200411073 

284.7312433 

30 

20011303 

3001014200 

SO 

9201IS 

13130739 

1331036 

322311337 

30 

920127 

270924SS 

1.731069 

230M6.73 

33 

0  5  10  IS  20  25  0  5  10  1  5  20  25 

3PEED  ANEMOMETER  in  m/s  SPEED  ANEMOMETER  in  m/s 

Rguic  5.8:  Mean  horiimntal  wind  speed  u  denved  wilfa  the  lidar  vemu  the  wind  speed  as  measured 
with  the  cap  anemometer  (kft).  Tlie  different  symbob  refer  to  different  data  files  as 
indicated  in  TaUe  5.4.  Ibe  mean  values  and  die  standard  deviatkais  (right).  Note  diat  the 
laigest  enor  ban  dominate. 

The  lelatioQ  between  die  borizontal  wind  diiectton  as  denved  widi  die  KAar  ud  naeasued  with 
the  wind  vane,  selected  from  different  data  files,  is  shown  in  Rgnte  S.9.  The  left  fignte  shows  die 
mean  values  of  different  data  files  and  different  periods;  the  right  fignte  shows  die  standard 
deviations. 
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«1A203*4n5 
•  e«7«8a9 


DIRECTION  ANEMOMETER  in  DEO  wrt  N 


♦ 


SPEED  ANEMOMETER  in  DEO  wrt  N 


ngDreS.9:  HorixoiMal  wnd  diractkn  as  deiivwd  wUi  tlw  IMv  Venn  the  wnd  direct 

wtt  llie  wind  vm  (left).  Difltapeal  qrnbols  lefler  to  dWenot  daM  files  as  indicated  in  TaUe 

^Same  as  the  right  bgl  with  sundaid  derisiioos  (left).  Note  that  the  iMgest  enor 

nan  tinanmBtM  " 


Horizontal  tneasotanaits  with  die  dual  Udar  have  shown  diat  in  12%  of  die  cases  dme  is 
snfiBdent  stroctore  for  the  deteiminatioo  of  die  wind  vector.  In  die  odier  cases  the  atmosphere 
was  weD  mixed  or  die  signal  to  noise  ratio  was  too  low.  A  «»"*»>  of  dual  KAmr  data 

files  widi  snfficieiit  stractme  were  analysed  for  cwirnUting  dw  wind  qieed  and  dhectioa.  b 

appeared  that  die  wind  dinctioa  as  derived  widi  die  lidar  dwwa  a  good  condation  with  the  data 

fiom  die  wind  vane.  On  die  odier  hand,  die  wind  speed  is  atarays  higbw  *«««  the  leadhigs  of  die 

ciq>  anemometer.  (This  was  a  reason  to  oompsae  the  cap  aaemometcrs  widi  Ml  odier  one.  b  a  later 

■tage.  die  cap  anemometer  was  compand  with  an  nltm  sonic  anemometer.)  As  expected,  in 

sitaatioos  when  the  wind  gpeed  oonld  also  be  dnnved  tram  the  slope  of  the  stnctnies  m  the  fidse 
oolonr  fignms.  the  renilts  txmespoito  wen  widi  diosB  obtaiaBd  with  dto  taoss^oiraiaiion  method. 
A  reason  fiv  the  discrepancies  between  the  lidar  and  the  cap  aicanomeicr  might  be  the  vertical 
component  in  die  wind.  Aaoendiag  and  deacendhig  stmctnres  can  dso  indnce  atmctnms  in  die 

horizontal  data  which  cannot  be  diatfaignished  from  the  horizontal  drift  of  the  stroctores.  This 
effect  has  been  indicaied  in  Rgnm  S.IO. 
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HORIZONTAL  ORirr  VERTICAL  DRIFT 


FfineS.lO:  Exinvle  of  M  abMonlMik  mcMe  wUch  ctmet  die  nme  dfect  in  a  Udar  dtta  let  if  it 
dtUb  eidier  horinatuvonreRkd^ 

An  odKr  reaicm  for  die  difcrapn(7  between  die  IkSv  nd  die  ciq>  aaemomBter  is  Ae  dififetc^ 
hei^  of  die  vobmiet.  Tbe  TNO  metecMower,  widi  die  cq>  anemometer  on  top,  is 

ffifuRM  in  a  hoDow  of  die  dunes  wbtle  die  lidar,  for  safety  reasons,  measures  ahntys  over  tbe 
dunes.  Because,  in  die  lowest  part  of  tbe  mixed  layer,  die  wind  incieases  as  a  function  of  altitude 
the  iviar  derived  wind  qieed  is  always  larger.  At  Ibis  stage  of  die  ejtycriment  it  is  not  possible  to 
compensate  for  diis  height  difference. 

5.4  Some  results  from  vertical  measurements 

Some  exanqdes  of  Bdv  returns,  measured  at  an  ekvation  angle  of  33  d^iees,  are  presented  in 
dus  section.  Bodi  die  taw  lidar  data  and  the  processed  data  are  shown,  including  die  actual  meteo. 
The  processed  data  is  obtained  from  the  raw  data  by  detienting  per  tow.  The  first  two  exaniples 
are  sdected  from  die  grotty  in  wfaidi  die  retnms  diow  onty  very  limited  amount  of  structure  at 
low  altitudes  while  die  odier  two  have  modi  more  structure  over  the  whole  vertical 

range.  Vertical  Bdar  measuremenls  provide  biformation  on  the  atmotyheric  structure  as  fimction 
of  altitude  being  die  basic  information  for  wind  calculatioos, 

Figure  5.11  an  exarafde  of  a  sitnarloo  where  die  atmotyhem  is  well  mixed  to  an  of 

dmot  700  m.  Above  this  hei^t  a  dood  is  detected,  uddcfa  provides  some  structure  in  die  interval 
from  700  m  to  900 m. 
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Pile:  911003-03122439-1224-1227 


MatBS: 

TempCTitnie  -  16,20  ±  0,02  *C 

Wetbalb  ±  *C 

ReLHmiidi^  -  77,09  ±  0,202  % 

PteMDie  -  1004,14  ±  7,136  liFi 
-  8,66  ±  0,19  km 


Tnmmiwkin  -  91,09  ±  1,21  % 

SdntiltatioD  -  11,19  ±  0,945  - 

^md^ieed  -  6,47  ±1,26  mrt 

WiDddfaectioii  -  224,36  ±  4,81  Degr. 


FigiaeS.12:  ReiidttoClitemewtHaBenlsttaBelevmlaoaqgleof 33dBgiees.lkiiintatkig8iiiiciavMaie 
detected  fram  dwot  300  to  Aoot  700  m  rittode.  Horinmial  azii:  0-187  ■;  veftlcal  azi«: 
0-1536111,  for  emAlkiv. 
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Some  moie  stnictuie  over  Ihe  i^Kde  range  of  the  lid«  is  shown  in  Rgnre  S.12.  However,  in  diis 
example,  noise  obscnres  die  atmoqibetic  stractnre. 

The  vrfiite  vertical  lines  in  the  lower  panels  of  Rgme  S.12  are  caused  by  imerfdence  from  die 
graphical  diqilay  to  die  Udar  receiver.  This  occurs  atmtimtitnedi  The  reason  for  this  cross-talk  is  not 
understood. 


S.S  Some  results  of  triangulatkin  measurements 

The  tr^le  lidar  has  been  »tniiiiar«>d  by  pointing  a  single  lidar  consecutively  in  three  different 
directions.  Hve  single  shot  measurements  in  eadi  direction  were  averaged  to  inqvove  the  signal- 
to-noise  ratio  and  to  tnwmniw  die  time  q)ent  to  point  die  lidar  in  die  desired  directioiis.  The  raw 
tr^k-lidar-data-set  of  die  first  exari^le  is  shown  in  Hgure  S.lSa,  the  data  after  detrendng  is 
shown  in  Figure  S.13b. 
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FignreS.lS: 

Rew  Ikter  date  ia  *«»«»>  ccdoor  m  nteMarod  te  three  different  direcikitt  in  a  ttere  venae 

heii^  ligiire  (lO-  Reealti  of  the  same  data  set  after  detrenting  (b).  Harizontal  axis:  0-1004  s; 
veitical  axis:  0u21<W  ol 
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The  data  set  shown  in  Rgme  S.13  was  analysed  wmi  the  assnnption  that  tibe  wind  vector  was 
stratified  horizontally.  No  vertical  conqxment  has  been  calcnlaled  here.  The  horiztwtal  wind 
direction  dins  derived  as  a  fhntrioo  of  altimde  is  shown  in  Rguie  S.14a.  The  wind  speed  as  a 
function  of  altitude  is  riiown  in  figure  5.14b.  The  variatkn  in  the  consecutive  data  points  can  be 
considered  as  an  indicatioa  for  die  uncertainty.  (In  dus  stage  of  data  processing  no  etnnr  bars  were 
calculated.)  A  second  exanqile  of  die  tr^de  lidar  method  is  presented  in  Rgure  5.  IS. 
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WIND  DIRECTION  vs  NORTH 


WIND  SPEED  in  m/s 


Hgrar  S.14:  Horizontal  wind  directkn  (left}  and  wind  qieed  as  a  functioQ  of  (right)  as  derived 

widi  a  triple  Udar  measaremem. 
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Mieteo: 

Teuipcratute 

-  3^1 

±  0,11 

•c 

Transmission 

-  90,98 

±  1,22 

% 

Wet  bulb 

-  2,01 

±  0,04 

•c 

Scmtillation 

-  6,29 

±  0,71 

- 

ReLHnmidily 

-  76,08 

±  1,59 

% 

Wind^eed 

-  0,87 

±  037 
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nessme 

-  1037,6 

±  0,21 

hpa 

Wind  direction 

-  280,86 

±  139 

Dbgr* 

Visibilior 

-  39,73 

±  0,77 

km 

Hgme  5.15:  Raw  H4<w  data  in  Mat  colour  as  measmed  in  tliiee  diCEoent  directions  in  a  time  versos 
heigiit  figure  (a).  Results  of  tiie  —ma  data  set  after  detrenting  (b).  Iforizootal  axis:  0-789  s; 
vertical  axis:  0-1440  m. 


TNOraport 


56 


The  apphed  algoridim  provided  no  error  bars  but  the  vatiadon  in  die  results  may  give  an 
impresskm  of  diis  aspect  At  altitudes  below  100  m,  die  results  were  calculated  at  eadi  available 
range  Un.  Above  100  m.  die  wind  vector  was  calculated  each  fifth  range  bin.  At  low  altitudes,  the 
calculated  wind  ditectioo  cone^mnds  widi  the  results  of  die  wind  vane.  At  higfMr  altitudes,  the 
qiread  in  die  wind  drrectkm  was  kas  than  at  lower  altitudes.  At  low  altitudes,  the  mlmlatarf  wind 
qmd  is  about  0.3  rrt/s.  This  is  widiin  die  range  of  die  wind  speed  as  measured  widi  the 
anemometer.  The  calailated  vertical  wind  speed  shows  mote  variation  at  low  altitudes  which 
might  be  caused  by  turbulence  or  by  the  uncertainty  in  the  data.  At  low  altitudes,  die  uncertainty 
in  («e  range  bin  or  one  time  bin  has  more  effect  than  at  larger  altitudes. 


No  reference  instrument  was  available  to  verify  die  results. 
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6  SUMMARY  OF  ACnVniES  AND  RESULTS  IN  THE  FIRST  PERIOD 

This  tepost  covers  part  of  die  activities  of  die  project  DWARSWIND',  A90K696,  is  die  period 
JttUMiy  1991  to  Match  1992.  In  that  period,  a  large  amonnt  of  softwaie  has  been  deveh^ied  for 
die  new  hardware  (controller  for  the  eqiupmeat,  graphics  di^Iay  and  die  fast  Fourkr  board. 
SeeKunz,  1992a.).  RndiennaK  the  bleachaUe  dye  Q-switcb  of  die  laser  was  rqdaced  by  a 
normal  Pockels-ceL  This  resulted  in  better  laser  stabili^.  The  interference  filters  in  die  receives 
were  rqilaoed  by  types  with  a  smaller  bandwiddi  and  die  noise  performance  of  die  receivers  was 
analysed  and  inqvoved.  It  is  eiqiected  diat  die  electronic  noise  can  still  be  reduced  by  about  a 
fector  5  if  other  conqiooents  are  used.  (See  Kunz  and  Moeiman,  1992.)  In  die  second  half  of  1991 
the  hard-  and  software  were  sufBciently  reliable  to  start  measurements  on  a  routine  basis  for  diis 
project. 

It  ^ipeaied  that  only  part  of  die  measurements  showed  (sufficient) .  ^mctuie  to  calculate  die  wind 
vector.  The  amount  of  structure  has  been  determined  subjectively  irom  false  colour  images  of  the 
data  sets  (for  diis  purpose,  only  the  variations  within  diree  times  die  standard  deviation  per  range 
bin  are  displayed).  The  results  of  this  classification  were  presented  in  section  5.2. 

In  the  second  semester  of  1S191,  software  was  developed  to  catoilate  die  horizontal  wind  vector 
from  the  (horizontal)  dual  lidar  measurements  and  to  calculate  the  wind  vector  as  a  functum  of 
altitude  frmn  lidar  measurements  in  dnee  different  (elevated)  directions. 

The  results  from  die  horizontal  measurements  show  that  there  is  a  good  correlation  between  die 
lidar  derived  wind  directiao  and  die  wind  direcdoD  as  measured  with  a  wind  vane.  The  wind 
qieed  as  derived  widi  die  lidar  is  higher  than  the  results  fiom  die  in  situ  measurements.  Possible 
reasons  for  diis  discrqiaiicy  could  be  die  vertical  ctmqxment  of  the  drift  of  die  structures  as  have 
been  discussed  in  section  5.3,  are  die  position  of  die  cup  MMaimnMttw  in  the  dunes  and  the 
difference  in  altitude  of  die  sensors.  The  time  required  to  perform  a  i«ti»r  measurement  is  of  the 
order  of  a  minute.  However,  ralniilarmg  of  dje  horizontal  wind  vector,  fif**  several  minntes  per 
range  bin  at  this  moment. 

By  measuring  simultaneousty  in  dnee  different  (elevated)  directioiu,  die  wind  vector  can  be 
determined  as  a  function  of  altitude.  The  princ^le  described  in  die  literature  is  based  on  die 
correlation  of  time  series  of  die  backacattered  radiation  from  dnee  different  time  series  at  a  fixed 
height.  The  proposed  mediod,  however,  derives  die  ditee-dimensional  wind  vector  as  a  function 
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of  altitude.  The  tr^le  hht  mediod  has  been  by  pointing  die  lidar  subsequently  in  dme 

different  directions.  Two  sets  of  lidar  data  with  sufBcfent  stmctme  were  analysed;  the  first  <»e 
widi  die  assunqption  of  a  horixnitally  Gratified  wind  vector  and  die  second  one  widiout  diis 
limitation.  Aldioagh  the  results  at  low  altitude  are  oonparable  widi  the  data  from  the  cap 
anemometer/wind  vane,  a  direct  in  situ  comparistm  is  not  possdde  at  die  momrat  The  time 
requited  to  measure  a  tr^  set  of  lidar  data  varies  from  IS  minntes  to  about  30  minutes.  This 
time  can  be  inqvoved  1^  using  die  seamd  Ikfar  with  a  better  logaridanic  amplifier  or  by 
omstructing  a  lidar  with  three  axes.  At  diis  mnneot,  conelation  calculations  for  the  three- 
dimensional  wind  vector  as  a  function  of  altitude  requires  several  hours  of  conqioter  time. 

It  is  worthwhile  to  verify  whether  lidar  ngnals  revealing  structures  are  (partly)  due  to  refractive 
index  variations  of  the  air.  A  brief  preliminary  discussion  of  diis  topic  is  presented  in  .^ipendix 
D.  In  die  framework  of  diis  project,  the  audior  has  described  this  topic  in  more  detail  in  a  sqiarate 
rqxnt  (Kunz,  1992c).  See  also  Kunz,  1993. 
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7  FUTURE  ACnvmES 

The  experimeiiis  peifofined  «id  the  lesulhi  cAtamed  have  shown  diat  die  wind  vector  can  be 
sensed  lemotely  widi  a  Mie  lidar  using  direct  detection.  The  technique  is  mainly  based  on  cross 
correlating  time  series  measived  at  different  positions.  Among  odiers.  a  task  for  the  future  is  to 
inqaove  dw  inversion  schemes  and  investigating  other  mathrmatical  tools  like  simulated 
Mingling  and  cross  spectral  phase  anatysis  as  proposed  by  Hooper  and  Eloranta  (1986). 

At  this  moment,  the  duee-point  measoietnent  is  based  on  measuring  consecutively  in  diiee 
different  directions.  With  die  dual  lidar,  however,  measurements  in  only  two  different  pointing 
directions  of  die  platform  will  sufBce  in  dieory. 

The  lidar  hardware  itself  can  also  be  improved  on  a  nundier  of  points: 

suppression  of  the  interference  udiich  occurs  during  die  rotaticHi  of  the  platform 

Stabilisation  of  die  laser  ouqmt  widi  a  new  Pockels-cel 

using  a  better  log-anqilifier  in  die  auxiliary  lidar 

inqirovement  of  the  mechanical  stability  of  the  beam  of  the  auxiliary  lidar. 

In  the  near  future,  lidar  measurements  wUl  be  performed  near  die  213  m  high  mast  of  the  Royal 
Meteorological  Institute. 
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HORIZONTAL  AND  VERTICAL  TANGENTS  TO  AN  OBUQUE  ELLIPSE 

This  qipendix  describes  some  bosks  for  snolysing  the  characteristk  ellipse  as  shown  in  Appendix 
B  and  C.  ifere,  only  die  coordinates  of  die  horixontal  and  vertical  tangents  to  an  oblique  ellipse 
are  calculated.  The  geometty  is  shown  in  Hgnre  A.1. 


Figure  Al:  Oblique  dlipse  with  its  horizontal  sad  votical  tangents. 

The  oblique  ellqise  diown  in  Hgure  A.1  can  be  described  by: 

(x  +  Vy)2  y2 


(A.1) 


wbett: 

2a  IB  lengdi  of  die  horizoiital  axis 
2b  «  kagdi  of  dw  vertkal  axis 
V  B>  coefficient  of  rotatka  of  the  eOipBe 

The  (iqiper)  horiztnital  tangoit  and  die  (right)  vertical  tangent  to  the  ellqise  are  determined  by 
consideting  equatwo  (A.  1)  as  a  inqilkU  function  of  x  and  y  which  can  be  described  as: 
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F(x,y): 


(x  +  Vy)2  y2 

- = -  +  -r  -  1 

*2  b2 


The  horizootal  tangent  follows  from: 


3F 

dx 


0 


which  leads  to: 

X  +  V  •y  =  0 

Substitution  of  (A1.4)  in  (Al.l)  gives: 

Xnj,  =  ii>V  and  yn*  =  ±b 
The  vertical  tangent  follows  from: 


9F 


(A.2) 


(A.3) 


(A.4) 


(A.5) 


(A.8) 


which  leads,  after  scmie  algebra,  to  the  line  on  which  the  extremum  must  lie: 


-Vb2x 
a2  +  v2b2 


(A.7) 


Substitution  of  (A.7)  in  (A.1)  leads,  after  some  algebra,  to: 

x^  »  ±V(a2  +  v2  b2)  (A.8a) 

and 


Vb2 

'^  +  V2b2) 


ymv  “  ± 


(A.8b) 
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DERIVATION  OF  THE  CHARACTERISTIC  PARAMETERS  OF  AN  OBUQUE  ELLIPSE 

An  oblique  ellipse  can  be  described  by  equadon  (B.l).  In  diis  q)pendix,  some  methods  aie 
described  to  calculate  die  parameters  a,  b  and  V. 

Apprwgh  1 


Input  parameters  are  the  intercept  with  the  vertical  axis  at  coordinates  (O.yg)  and  the  vertical 
tangent  at  cowdinates  (xi,yj). 


Hgure  B.l:  An  oblique  ellipse  with  the  intercept  at  coordinates  at  (0,y(^  and  the  vertical  tangent  at 
coordinates  (X},y|). 


The  oblique  ellqise  is  described  by: 


(x  +  Vy)2  y2 

- = - +-=•  *  1 

t?  b^ 


(B.l) 


where: 

2a  s  length  of  the  horiztmtal  axis 
2b  s  length  of  the  vertical  axis 
V  a  coefficient  of  rotation  of  the  ellqise 


The  ellipse  is  shown  with  the  horizontal  and  vertical  tangents  in  Figure  B.l. 
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The  basic  equations  to  calculate  the  parameten  a,  b  and  V,  of  tite  ellqwe  are: 


coordinate  (0,yo) 

(B.2) 

Xi2  =  a2  +  b2  ■  V2,  see  equation  (A.8a) 

(B.3) 

yi  -V 

Xi  “  (a2  +  v2  b2)’ 
Solving  parameter  b: 


see  equation  (A.8b) 


(B.4) 


From  (B.3)  and  (B.4)  an  expiessitm  for  ■  b^  is  derived.  This  leads  to: 


V2  1)2  » 


yi2xi2 


(B3) 


Using  (B.l),  (B.2)  and  (B.3)  an  other  ejqnessicm  for  V2  •  b2  is  derived: 

Xi2  (b2.yo2) 


V2  .b2 


b2 


(B.6) 


Combining  the  equatirms  (B.S)  and  (B.6)  leads  to: 

b  =  V(yo2  +  yi2) 

Solving  die  parameter  V: 

Substitution  of  (B.7)  in  (B.5)  provides  an  exptessim  for  V: 

y-  yi-^1 

(yo*+yi*) 

Solving  fee  pfffaniptfir  a; 


(B.7) 


(B.8) 


Substitution  of  (B.7)  and  (B.8)  in  (B.3)  provides  an  eqnession  for  a: 


TNOrapoit 


AppcBdix  B 


P«gc 

B.3 


a  =  *iyo 
^(yi^  +  yo^) 


(B.9) 


rnncliMinn: 

Hie  parameters  of  the  characteristic  ell^ise  can  be  calculated,  if  the  coordinates  of  die  intercept 
with  die  vertical  axis  and  the  coordinates  of  the  vertical  tangent  are  given. 

Approach! 

The  characteristic  parameters  of  an  oblique  ellqise  can  also  be  determined  if  the  intercepts  with 
the  vertical  axis  and  widi  the  horizontal  axis  are  given  and  die  intercqit  widi  die  vertical  dirough 
the  horizontal  intercept.  This  is  shown  in  Hgure  B.2. 


Figure  B.2:  The  chaiactaistic  ellipse,  croning  the  vertical  and  the  horizontal  axis  renectivefy  at 
®o<*<bnates  (0,y^  and  (xj,0).  A  thud  coordiiiate  is  given  by  the  vertical  line  dirough  (xn,0) 
adiidi  crosses  the  ellipse  at  (X|,y|).  ^ 


The  equations  to  calculate  die  parameters  of  die  oblique  ellqise  ate: 


y-0 


(B.IO) 
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(Xi  +  Vyi)2  yi^ 


(B.ll) 


v^  ya^  yjS 

a2  ■*■52 


1. 


(0,y2) 


(B.12) 


Substituting  of  (B.IO)  in  bodi  (B.ll)  and  (B.12),  while  eliminating  b2  results  in  an  e]q>licit 
equation  for  V: 


(Xi+Vy,)2  1  V2  1 


yi^  *1^  ya^ 


(B.13) 


Some  algebra  leads  to: 


(B.14) 


The  parameter  b  can  be  found  by  die  subsdtuticm  of  equations  (A2.10)  and  (B.14)  in  (B.12).  This 
results  in: 


yr  1  _  1 

4ya^^b2 


(B.15) 


which  leads  directly  to  a  solution  for  b2: 


(B.16) 


T 


•w 
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A  GENERAL  SOLUTION  OF  THE  CHARACTERISTIC  PARAMETERS  OF  AN  OBUQUE 
ELLIPSE 


In  Appendix  B,  a  sohitioa  was  given  to  derive  the  characteristic  parameters  of  an  oblique  ellipse 
in  two  qwdal  situatUMis.  hi  this  qipendix,  an  attempt  will  be  made  to  solve  die  characteristic 
parameters  in  die  case  diat  die  intercept  widi  die  vertical  axis  is  given  and  two  odier  coordinates, 
which  describe  die  intercqits  widi  an  arbitrary  vertical,  as  shown  in  Figure  C.l. 


Figure  Cl:  The  chaiactenstic  ellqise  widi  die  intercept  at  the  votical  axis  and  two  points  crotsii^  an 
arbilraiy  verticaL 

In  princ^le,  die  problem  can  be  solved  coixqiletefy  because  three  parameters  describe  the  ellipse 
and  dnee  coordinates  are  ghren.  Ihiwever,  it  appeared  diat  dus  problem  can  only  be  solved  widi  a 
very  lengthy  algdna  whidi  has  not  fully  been  described  here.  For  example,  die  final  equation  for 
die  parameter  V  is: 


I 


1 
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-2  xi-  y2  2  xi  2  X|  2  xi 

ya  ya*  yrya^  yrya^ 


2 

yi^  ya^  ya^  ya^ 

The  solutions  for  the  parameters  b  and  a  are  even  more  coiiq>]icated. 

Although  diis  solution  is  more  general  and  includes  the  odier  two  described  in  Appendix  B,  the 
con^lexity  of  the  solution  does  not  justify  its  elidxnation  at  this  moment. 

Nevertheless,  dte  advantage  of  diis  method  would  be  that  a  more  general  solution  becomes 
available  which  can  use  differmt  values  of  the  cross-  and  autocorrelation  functions  as  discussed 
in  Chapter  4  and  5  of  diis  report.  This  would  inprove  the  quality  of  the  solution  because  a  larger 
interval  of  the  functions  is  used  in  stead  of  only  one  point. 
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AN  ESTIMATION  OF  REFLECnCW  OF  AIR  LA  YERS  WITH  DIFFERENT  TEMPERATURE 
AND  PRESSURE 

Thus  te,  it  has  been  in  dus  lepoct  (and  also  in  the  literature)  that  atmo^theric  structures, 

as  measured  widi  lidar,  are  a  result  of  (drifting)  aetoatd  structures.  The  differeiice  in  aerosol  type 
and/or  concentratiaii  of  structures  causes  a  modulation  (m  the  badcscitttered  signal. 
However,  it  is  also  posyiWe  diat  the  signal  varies  due  to  local  differences  in  refractive  index.  The 
amooqiheric  refractive  index  varies  with  tenqterature,  pressure  and  relative  humidity.  If  diese 
variations  are  large  «^«ongh,  di^  can  act  as  i^Mtial  Iresnel  reflectors. 

The  objective  of  diis  section  is  to  investigate  the  effects  of  local  refractive  index  variations  tm  the 
reflected  signal. 

The  value  of  the  refractive  index  of  air,  wfafeh  varies  with  tenqxrature  and  pressure,  can  be  found 
in  several  publicatitms  like  Btooknei(1977),  Cliffott  (1977),  Zuev  (1982)  and  Falcone  and  Dyer 
(1985). 


(DL-  1)  =  77.6  (1  +7.52  •  10-3  .  X-2)  .  10^ 

T 


(D.l) 


in  which; 

n,  >  refractive  index  of  air 
X  m  waveiengdi  of  die  radiation  in  pm 
P  s  pressure  in  millibars 
T  >  temperature  in  K 

Literature  on  dus  subject  indicates  diat  the  relative  humidity  has  only  a  minor  infhience  on  the 
refractive  index  for  radiatkm  in  die  visible  part  of  die  electromagnetic  spectrum.  Therefore,  die 
total  variation  in  die  refractive  index  can  be  described  by: 
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In  an  atmoq>hete  with  a  teirtprramre  of  300  K  and  a  paeaaise  of  1000  mBar,  the  refractive  index 
at  a  wavekngdi  of  1  fun  variea  aboot  10*^  if  the  tempenttae  varies  0.1  K  and  die  pressure  varies 
0.1  mBar.  As  a  result,  the  reflectioa  is  about  lO'K  hr  our  situation,  where  the  laser  power  is  aboot 
1  MW  and  the  noise  equivakat  power  of  die  receiver  is  about  5  nW,  die  amount  of  power 
reflected  from  such  an  air  layer  is  detectable. 

CfliKhmiffli 

A  first  investigatkn  shows  diat  step  like  difietences  in  temperature  and  pressure  might  introduce 
variations  in  refractive  index  of  air  whidi  are  strong  enough  to  produce  detectable  reflections  in 
the  lidar  receiver.  A  trK»e  extended  description  of  diis  subject  has  been  published  by  die  audior 
(Kunz,  1992c). 


ONGERUBRtCEERD 


REPORT  DOCUMENTATION  PAGE  (MOO-NL) 


1.  DEFENSE  REP0m'NUMBfR(M004IL)  2.  RECIPiBirS  ACCESSION  NUyBER  S.  PBVOWmWIQ  ORGANIZATION  REPORT 

NUMBER 

TD93-0485  FEL-93-A040 


4.  PflOJECT/TASK/WORK  UMT  NO. 
22393 


5.  CONTRACT  NUMBER 
A90K696 


E  REPORT  DATE 
AUGUST  1993 


7.  NUMBER  OF  PAGES  8.  NUMBER  OF  REFERENCES 

75  (INCL  4  APPENDICES,  39 

EXCL  RDP  +  DISTRIBUTION  LIST) 

t 


9.  TYPE  OF  REPORT  AND  OATES  COVERED 


10.  TITLE  AND  SUBTITLE 

WIND  MEASUREMENT  WITH  AN  INCOHERENT  UDAR  (RRST  PROGRESS  REPORT) 


11.AUTHOR(S) 
G.J.  KUNZ 


PERPORMMO  ORGANIZATION  NAME(S)  AND  AOORESS<ES) 

TNO  PHYSICS  AND  ELECTRONICS  LABORATORY,  P.O.  BOX  96864, 2509  JG  THE  HAGUE 
OUDE  WAALSDORPERWEG  63,  THE  HAGUE,  THE  NETHERLANT^ 


13.  SPONSORINGAMNITORING  AGENCY  NAME(S) 
ROYAL  NETHERLANDS  NAVY 


14.  SUPPLB4ENTARY  NOTES 

THE  CLASSIFICATION  DESIGNATION  ONGERUBRICEERD  IS  EQUIVALENT  TO  UNCLASSIFIED. 


IS.  ABSTRACT  (MAXIMUM  200  WORDS,  1044  POSITIONS} 

ATMOSPHERIC  STRUCTURES  CAN  BE  MEASURED  WITH  INCOHERENT  OPTICAL  RADARS  (UDARS).  BECAUSE  THESE 
STRUCTURES  DRIR  WITH  THE  WIND,  THEY  CAN  SERVE  AS  A  TRACER  FOR  REMOTE  SENSING  OF  THE  WIND  VECTOR.  FOR 
THIS  PURPOSE,  A  DUAL  MONOSTATIC  SCANNING  UDAR  SYSTEM  IS  AVAILABLE  TO  MEASURE  THE  ATMOSPHERE 
SIMULTANEOUSLY  IN  TWO  DIFFERENT  DIRECTIONS  OVER  A  MAXIMUM  RANGE  OF  ABOUT  1  KM.  THE  TRANSIT  TIME  OF 
IDENTIFIED  PATTERNS  BETWEEN  TWO  SENSING  POINTS  IN  THE  HORIZONTAL  PLANE  PROVIDES  IN  COMBINATION  WITH 
THE  GEOMETRY  OF  THE  UDAR,  SUFFICIENT  INFORMATION  TO  DERIVE  THE  HORIZONTAL  WIND  VECTOR.  THE  METHOD  IS 
BASED  ON  CROSS-CORRELATION  TECHNIQUES.  TO  DETERMINE  THE  SPATIAL  WIND  VECTOR  AS  A  FUNCTION  OF 
ALTITUDE  IT  IS  SUFFICIENT  TO  MEASURE  IN  THREE  DIFFERENT  UPWARD  DIRECTIONS.  THIS  CAN  BE  REALIZED  WITH  A 
TRIPLE  UDAR  OR  WITH  A  SINGLE  UDAR  BY  MEASURING  CONSECUTIVELY  IN  THREE  DIFFERENT  DIRECTIONS  AND  USING 
AN  EQUIVALENT  BUT  MORE  EXTENDED  INVERSION  METHOD.  THIS  REPORT  DESCRIBES  THE  MAIN  ACTIVITIES  CARRIED 
OUT  FOR  THE  PROJECT  DWARSWIND’  IN  THE  PERIOD  JANUARY  1991  -  MARCH  1992.  AFTER  A  SUMMARY  OF  A  SHORT 
LITERATURE  STUDY,  SOME  THEORETICAL  ASPECTS  ARE  DESCRIBED  SUCH  AS:  THE  VECTOR  REPRESENTATION  OF  THE 
DUAL  UDAR  IN  CARTESIAN  COORDINATES,  A  METHOD  TO  DERIVE  THE  WIND  VECTOR,  THE  CHARACTERISTIC  LIFE  TIME 
AND  THE  CHARACTERISTIC  SIZE  OF  THE  STRUCTURES  FROM  A  SET  OF  LDAR  MEASUREMENTS.  A  SELECTION  OF 
EXPERIMENTAL  REaiLTS  ARE  PRESENTED.  FOR  EXAMPLE,  THE  CALCULATION  OF  THE  HORIZONTAL  WIND  VECTOR  AT 
AN  ALTinOE  OF  15  M  AND  THE  WIND  VECTOR  UP  TO  AN  ALTITUDE  OF  1000  M. 


IS.  DESCmPTORS  IDBmRERS 

WIND  MEASUREMENTS 
LASER  RADAR  (IX>AR) 

VERTICAL  DISTRIBUTION  OF  WIND 
HORIZONTAL  DtSTRSUTION  OF  WIND 
WIND  VELOCITY  METER  (ANEMOMETER) 


17a  SECURITY  CtASSIFICATION 
(OP  REPORT) 
ONGERUBMCEERD 


17bk  SECURITY  CLASSIFICATION 
(OF  PAGE) 
ONGERUBRICEERD 


17a  SECURITY  CLASSIFICATION 

(OF  ABSTRACT) 
ONGERUBRICEERD 


IS,  DISTRIBUTIONIAVAILABSJTY  8TATEMBIT 
UNLIMITED 


17d.  SECURITY  CLASSIFICATION 
(OF  TITLES) 
ONGERUBRICEERD 


ONGERUBRICEERD 


